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ABSTRACT 

We present a theoretical model for Type lb supernova (SN) 2006jc. We calculate the evolution of the 
progenitor star, hydrodynamics and nucleosynthesis of the SN explosion, and the SN bolometric light 
curve (LC). The synthetic bolometric LC is compared with the observed bolometric LC constructed 
by integrating the UV, optical, near- infrared (NIR), and mid- infrared (MIR) fluxes. The progenitor 
is assumed to be as massive as 40Mq on the zero-age main-sequence. The star undergoes extensive 
mass loss to reduce its mass down to as small as 6.9Mq, thus becoming a WCO Wolf-Rayet star. The 
WCO star model has a thick carbon-rich layer, in which amorphous carbon grains can be formed. This 
could explain the NIR brightening and the dust feature seen in the MIR spectrum. We suggest that 
the progenitor of SN 2006jc is a WCO Wolf-Rayet star having undergone strong mass loss and such 
massive stars are the important sites of dust formation. We derive the parameters of the explosion 
model in order to reproduce the bolometric LC of SN 2006jc by the radioactive decays: the ejecta 
mass 4.9M0, hypernova-like explosion energy 10^^ ergs, and ejected ^^Ni mass O.22M0. We also 
calculate the circumstellar interaction and find that a CSM with a flat density structure is required 
to reproduce the X-ray LC of SN 2006jc. This suggests a drastic change of the mass-loss rate and/or 
the wind velocity that is consistent with the past luminous blue variable (LBV)-like event. 
Subject headings: dust, extinction — infrared: ISM — nuclear reactions, nucleosynthesis, abundances 
— supernovae: general — supernovae: individual (SN 2006jc) — stars: Wolf-Rayet 



1. INTRODUCTION 

On 9th October 2006, iNakano et all (|2006[ ) reported 
K. Itagaki's discovery of a possible supernova (SN) in 
UGC 4904. Although the SN was discovered after 
the peak, an upper limit of the R magnitude (M/j > 
— 12 .2) was obtained at ~20 days befor e the discov- 
ery (iPastorello et al.ll2007f ). Interestinglv. iNakano et all 
()2006( l also reported that an optical transient had ap- 
peared in 2004 close to the position of SN 2006jc. The 
transient was as faint as Mr ^ — 14 and its duration 
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was as short as ^ 10 days. Since the event was faint 
and short-lived, they speculated that the transient was 
a luminous blue variable (LBV)-like event. The spatial 
coincidence bet ween the LBV-like eve nt and SN 2006jc 
is confirmed bv IPastorello et aP ()2007| ). Because of such 
an intriguing association with the LBV-like event, many 
groups performed follow-up observations of SN 2006jc in 
various wavebands: X-ray, ultra violet (UV), optical, in- 
frared (IR), and radio. 

Spectroscopic observations showed many broad fea- 
tures and strong narrow He I emission lines. According 
to the He detection, S N 2006.ic was classifie d as Type 
lb (ICrotts et al. 20 06: F esen et a l. 2006a; Be netti et al.l 
l2006t iModiaz et aLi 2006a, b). However, strange spectral 
features and their evolutions were reported. A bright 
blue continuum was promiiient in the optical spectruni 
at early epochs (Fol ev et al]|2007t IPastorello et al.ll2007t 
iSmith et a l. 2008). Such a bright blue c ontinuum had 
also b een observed in Type II SN 1988Z (jTuratto et all 
Il993| ). but the origin of this feature is still unclear. As the 
blue continuum declined, the red wing bri ghtened and 
the optical spe ctra showed "U"-like shapes ("Smith et al.l 
2008; Kawabat a et al.ll2008D . This is a distinguishing fea- 
ture of SN 2006JC in contrast to the spectra of usual SNe 
that have a peak in optical bands. 

Photometric observations in optical and IR bands were 
performed continuously. The optical light curve (LC) 
showed a rapid decline from 50 days afte r the disco very, 
as in the case of SN 1999cq (jMatheson et al.ll2000h . At 
the s ame epoch, near infra r ed (NIR) emissions bright- 
ened (jArkharov et al.|[200l iDi Carlo et al.ll2007D . The 
NIR brightness increased from ~ 40 day s to ^ 70 days 
after the discovery and then declined (jDi Carlo et al.l 
|2007| ). The epoch of the NIR brightening corresponds 
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TABLE 1 
Optical luminosities. 
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Fig. 1. — Multicolor available observations of SN 2006jc. (a) 
X-ray and UV luminosities obtained with the Swift and Chan- 
dra satellites (black: X- ray, red: UV W2-band, green: UVM2- 
band, blue: UVWl-band. Ilmmler et al . 2008). X-ray LC is shown 
in a unit of ergs s~^ (right axis), (b) Optical luminosities ob- 
tained with the HCT telescop e (fille d, Anupama et al. 2008) and 
the SLTBARU" telescope (open, Kawabata ct al. 2007al. 120081) . The 
upper limit is taken from j'astorcllo et al. (200'5). The color rep- 
resents the wavelengths (U-band: red, B-band: green, V-band: 
blue, R-band: magenta, I- band: cyan), (c) IR luminosities ob- 
tain ed with the AZT-24 telescope (filled circles, lArkharov et al.l 
[2006: Di Carlo ct al. 2007), the MAGNUM telescope (open cir- 
cles, Minczaki c t al.i i2007D . and the AKARI satellite (squares, 
ISakon ct al. 2007). The color of circles represents the wavelengths 
(J-band: red, H-band: green, and K-band: blue). The contribu- 
tions to the IR luminosities from the hot dust (filled square) and 
with the hot and warm dust (open square) are shown in a unit of 
ergs s~^ (right axis). 

to that of the development of the red wing in the optical 
spectra (jSmithet al.ll2008f) . 

The NIR brightening, as well as the fact that the redder 
side of the He emission profile declined faster than the 
bluer side, has been interpreted as an evidence of an on- 
going dust formation (Smith ct al. 2008). Additionally, 
on 29th April 2007 (200 days after the discovery), the 
AKARI satellite performed NIR and mid-in frared (MIR) 
photo metric and spectroscopic observations (jSakon et al.l 
120070 and the MAGNUM telescope obtained the NIR 
photometries (Minczaki ct al. 2007). They report the 
formation of amorphous carbon dust: another piece of 
evidences of the dust formation. 

X-ray and UV emissions have als o been observed by 
the Swift and Ch andra satellites (Brown e t al.l 120061 : 
Ilmmler et aLll2006l 2008; Holland ct al. 2007). X-ray ob- 
servations were performed at seven epochs and showed 
a brighten ing from ~ 20 days to ^ 100 days after th e 
discovery (jBrown et al.ll2006l: HiWler et al.ll20'oll2008D . 
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The X-ray detection suggests an interaction between the 
SN ejecta and the circumstellar matter (CSM). On the 
contrary, the radio emission was no t detected by Very 
Large Array (VLA) (|Soderberd l2006'). 

We present a SN explosion model of a Wolf-Rayet star 
that explains the bolometric and X-ray LCs. Hydrody- 
namics, nucleosynthesis, and LC synthesis calculations 
are performed assuming the spherical symmetry. In this 
study, we assume the explosion date of SN 2006jc to be 15 
days before the discovery (t = 0) and the energy source 
of the light to be the ^^Ni-^SCo decay. 

The paper is organized as follows: in § |2l we describe 
how we derive the bolometric LC from observations in 
the various wavebands, in §|31 we briefly discuss the pre- 
supernova evolutionary properties of the progenitor star; 
in § m hydrodynamical and nucleosynthesis calculations 
are described; in § |5l LC synthesis calculations are pre- 
sented; in § [6l we calculate the X-ray emission due to the 
ejecta-CSM interaction; in § |7| and § El conclusions and 
discussion are presented. 

2. PHOTOMETRIC OBSERVATIONS AND BOLOMETRIC 
LIGHT CURVE 

The bolometric luminosities of SNe are usually esti- 
mated from the integration over the optical and NIR 
emission because the usual SNe rad iate dominantly in 
the optical and NI R bands (e.g., lYoshii et alJ 120031 : 
iMinezaki et~ani2008l ) . However, the s pectra of SN 2006.ic 
show the bright red and blue wings (ISmith et al.|[2008l : 
iKawabata et al.ll2008l : I Anupama et al.ll2008l ). which im- 
plies that the emissions in UV and IR bands considerably 
contribute to the bolometric luminosity. 

We construct the bolometric luminosity with the inte- 
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Fig . 2 . — C omparison of multicolor LCs of SN 2006ic (Lppt: green 
filled circles, IKawabata et al.ll2007al. [2003: 1 Anuparaa et all l2008l: 
Luv + ^opt: cyan filled circles Mmmler et al .' 2008; Anupam a et al.l 
[2003: LiR,est(i^ < 3 X lO^^Hz): red filled circles. Arkharov et al.l 
[200^: IDi Carlo et"all \20m\. Lppt + i^iR. e st (t^ < 3 X lO^^Hz): 
blue filled circles. [ K awabat a et a l.l |2007a|, I200 8|: l Anupama et al.l 
[2003: lArkharov et al.l 120061 : IDi Carlo ct al. 2003; ilR,hPt(;^ < 
3 X lO^^Hz): magenta filled circle, .Sakon ct all |2007| : Li^{u < 
3xlO^*Hz): magenta open cire/eJSakon et al.l2007l : and Lx: black 
filled cirefe^. Flmmler et al.|[2003) . 
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scopes (IKawa bata ct al. 2007a, 120081 ) an d the Swift 
(llmniler et al.l 12008) and AKARI satellites ()Sakon et all 
|2007( ). Since the UV fluxe s are available only at i = 17 
days (jlmmler et al.l I2008D . the UV luminosity is esti- 
mated from the optical luminosity at the other epoch. 
Available observations are shown in Figure [TJ Details of 
optical observati ons will be presented in the forthcom- 
ing p apers (e.g., lAnupama et al.l [20081 : IKawabata et al.l 
|2008[) . We adopt a distance o f 25.8Mpc correspondi ng to 
a distance modulus of 32.05 (jPastore llo ct al. 2007^ and 
a reddening of E(B - V) = 0.05 dSclilcgcl ct al,, 1.998,; 
iPastorello et al.ll2007l ). 

2.1. Optical emission 

The optical LCs wer e obtained with the HCT 
and SUBARU telesco pes (jKawabata et al.ll2007aL 120081 : 
lAnupama et al.|[2008( ). We integrate the optical fluxes 
with a cubic spline interpolation from 3 x 10^^ Hz to 
1 X 10^^ Hz. The optical luminosities (iopt) are sum- 
marized in Table [1] and the LC is shown in Figure [3 
The optical LC declines monotonically after the discov- 
ery. The decline suddenly becomes rapid at i > 70 
days and the optical luminosity flnally goes down to 
Lopt ~ lO^^ergs s'^ at t ~ 200 days. 

The X- ray LC obtained w ith the Swift and Chandra 
satellites (jlmmler et al.l[2008f ) shows that the X-ray lu- 
minositie s, Lx, are much fa i nter than the optical lumi- 
nosities (jBrown et al.l 120061: llmmler et af] l2006l I2008D . 
Thus, the X-ray contribution to the bolometric lumi- 
nosities is negligible. However, the UV luminosity, iuv, 
is comparable to the optical luminosity at t = 17 days 
(Luv '^ 3 X lO'^^ergs s~^ as estimated from the UVOT 



observations. llmmler et aL[|2008l) .^^ The UV luminosity 
is ~ 80% of the optical luminosity, i.e., the total flux is 
~ 1.8 times brighter than the optical flux fFig.lH). Since 
the U V flux declined as the optical flux ( Holland et al.l 
12007( 1. we assume that iuv ^ 0.8Lopt at every epoch. Al- 
though the blue wing declines wi th time and Ltiv m ight 
be over-estimated at i > 90 days ([Smith et al.1120081 ) , the 
bolometric luminosity (Lboi) should be reliable because 
the IR contribution dominates in the bolometric lumi- 
nosity at such late epochs (iJ I2.2p . 



2.2. Infrared emission 

The IR spectroscopy and photometries are obtained 
with the A ZT-24 and MAGNUM telescopes (NIR pho- 
tometries, 'Arkhar ov et all [20061 : iDi Carlo et all [20071 : 
iMinezaki et al. 2007[) and the AK ARI satellite (NIR 
spectroscopy and MIR photometries. ISakon et al.ll2007f l. 
As indicated by the red wing in the optical spectra, the 
IR emission considerably contributes to the bolometric 
luminosity of SN 2006JC. 

The MIR observ ation is available at i = 215 days 
([Sakon et al.l [20071 ). The IR luminosity integrated over 
V <3x 10" Hz is estimated from the NIR and MIR ob- 
servations as Li^jv < 3 x lO^^'Hz) = 4.5 x lO''" ergs s^^. 
I Sakon et al.l ([20071 ) concluded that the IR emission is 
originated from amorphous carbon grains with two tem- 
peratures of T = 800K and 320K. The large difference 
between the two temperatures would imply that the ori- 
gin of the hot carbon dust with T = 800K is different 
from that of the warm carbon dust with T = 320K. 
The hot carbon dust is suggested to be newly formed 
in the SN ejecta and heated by the ^^Ni-^^Co d ecay by a 
dust formation calculation ([Nozawa et al.|[2008r ) . On the 
other hand, the origin of the emission from the warm car- 
bon dust is s uggested to b e a SN hg ht echo of the C SM 
car bon dust (ISakon et all [2007: Mat tila et al.l [20081: see 
also [Nozawa et al] |2008). Therefore, we assume that the 
optical emission from SN 2006jc is absorbed and simul- 
taneously re-emitted by the hot carbon dust and thus 
the luminosity emitted from the hot carbon dust should 
be included in the bolometric luminosity of SN 2006jc. 
According to the estimated temp eratures and mass es of 
the hot and warm carbon grains ([Sakon et al.ll2007[ ). the 
luminosities contributed by the hot and warm carbon 
grains are LiR^hot(i^ < 3 x lO^Hz) = 3.2 x 10"° ergs s"! 
and LiR,warm(i^ < 3 X lOi'^Hz) = 1.1 X 10"*" ergs s^^, 
respectively.^^ 

For the epochs when the IR photometries at z/ < 1.3 x 
lO^^'Hz are unavailable, we estimate the contribution of 
the IR emission by fitting the JHK-band photometries 
with amorphous carbon emission. 

From the Kirchhoff's law, the thermal radiation from 
a spherical dust grain X with a uniform radius ax and 
temperature Tx is given by ATTa^B{i', Tx)Qx^ (i^) , where 
Qx^i'^) is the absorption efficiency of the grain. For 
the optically thin case, the observed emission from dust 

^^ The UV flux is estimated using the Swift UVOT calibrat ion 
documents ( http://swift.gsfc.nasa.gov/docs/swift/swiftsaMmlb. 

^^ The difference between Li^(u < 3 X lO'^^Hz) and Ljr hptl'^ < 
3 X 10l-*Hz) + LiR,warm(!^ < 3 X lO^Hz) stems from that the H- 
band luminosity is slightl y brighter than th e luminosity emitted 
from the hot carbon dust llSakon et al.ll2007l ). 



TABLE 2 

Parameters for amorphous carbon fitting of the JHK-band photometries and 

the estimated ir luminosities. 
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grains X is written as 



fxii^) = Nx7rB{,y,Tx)Qi'' 






(1) 



where TVx and R denote tlie total number of the dust 
particles and the d istance from the observer, respectively 
(jSakon et al.ll2007| ). In the foUowings, we convolve the v- 
independent coefficients as an emission coefficient C^ = 
ttTVx (ax/-R) ■ Applying the absorption efficiency for the 
amorphous carbon grain with ac = 0.01/im, we derive 
the temperature of the hot carbon dust, Tchot, and C^ 
to reproduce the JHK-band photometries. 

To justify the above estimate, we compare the es- 
timate WTth__theactual MIR observation at t — 215 
days (jSakon et al.l l2007l ). The fitting gives the tem- 
perature Tc,hot = 870 K and the emission coefficient 
Ce = 2.8 X lO^^ for the HK-band photometries at t = 215 
days. The luminosity integrated over i^ < 3 x lO^'^Hz is 
U^,^Av < 3 X lOi-^Hz) = 2.7 X 10-*" ergs s"!. The 
temperature and luminosity are roughly consistent with 
those of the hot carbon dust. The agreement indicates 
that the fitting gives a good estimate of the IR emission 
due to the hot carbon dust. We note that the estimate 
can not account for the emission from the warm carbon 
dust. 

Table [2] summarizes the emission coefficient, temper- 
ature, estimated luminosity at ;/ < 1.3 x 10^'* Hz, and 
luminosity emitted below i/ = 3 x 10^^ Hz. The dust tem- 
perature roughly decfines from Tc,hot ~ 1600K at i = 49 
days to Tchot ~ 870K at t = 215 days. This is consistent 
with a picture that the hot carbon dust was formed in 
the S N ejecta and cooled down gradually (jNozawa et al.l 
|2008[ ). The IR LC is shown in Figure H The estimated 
luminosity at ly < 1.3 x 10^^ Hz evolves as the JHK LCs, 
and thus the IR LC brightens at i ~ 50 — 80 days and 
declines at i > 120 days. Since there is no NIR data 



TABLE 3 

BOLOMETRIC LUMINOSITIES. 
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at i ~ 80 — 120 days, the bolometric LC can not be 
estimated at this epoch. The bolometric luminosity is 
derived from the summation of iuVi -^opt, and Ljr and 
summarized in Table [3l where Luv = O.SLopt is applied. 

3. THE PROGENITOR STAR 

The presupernova model has been extracted from a set 
of models already presented bv lLimongi fc Chieffil (|2006D 
and computed with the latest release of the stellar evolu- 
tionary code FRANEC (5.050218). Since all the features 
of this code have been already presented, we will address 
here only the main points. The interaction between con- 
vection and local nuclear burning has been taken into 
account by coupling together and solving simultaneously 
the set of equations governing the chemical evolution due 
to the nuclear reactions and those describing the convec- 
tive mixing. More specifically, the convective mixing has 
been treated by means of a diffusion equation where the 
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Fig. 3. — Abundance distributions of presupernova models with (a) Afnis = 25A/q, (b) A/ms = SOA/q, (c) Afms = 4OAf0, and (d) 
Afms = 120AfQ. Note that the chemical composition of the outmost layer is C rich in the WCO Wolf-Rayet stars {c,d). 



diffusion coefficient is computed by the use of the mixing- 
length theo ry. The nuclear network: i s the same as that 
adopted in iLimongi fc Chieffil ()2003l ). but the nuclear 
cross secti ons have been updated w henever possible (see 
Table 1 in ILimongi fc Chieffill2006[ ). A moderate amount 
of overshooting of 0.2 Hp has been included into the cal- 
culation only on the top of the convective core during 
core H burning. Mass loss h as been taken into account 
following the prescriptions of IVink et alj (l2000f) for the 
blue s upcrgiant phase (Tes > 12000K), Ide Jager et al.1 
(119881) for the r ed supergiant phase (Tcs < 12000K), 
iNugis fc Lai^ (I2OOQI for the WNL Wolf-Rayet phase 
and ILange^ 119891) during the WNE/WCO Wolf-Rayet 



phases. We adopt the following correspondence of the 
models to the various WR pha ses according to the surfac e 
abundances, as suggested bv iMaeder fc Mevnetl (|2003D : 
WNL (10-5 < x(H),urf < 0.4), WNE (X(H)surf < 10"^ 
and (C/N),urf < 0.1), WNC (0.1 < (C/N)surf < 10) and 
WCO ((C/N)surf > 10). (Hereafter, C/N and C/0 de- 
note the number ratios.) 

The X-ray emission, as well as the early bright blue 
continuum and the narrow He I lines, clearly indicates an 
interaction between the SN ejecta and the CSM, i.e., the 
existence of a dense CSM. Furthermore, the IR spectral 
energy distribution may be explained by the formation 
of amorphous carbon grains in th e SN ejecta and th e 
CSM (Sakonetal. 2007i, see also iNozawa et all 120081 ). 
Since the C-rich environme nt (i.e., C/0 > 1) is required 
to form carbon dust (e.g., iNozawa et al.l 120031 ) . the IR 



observations suggest that the SN ejecta and CSM contain 
a C-rich layer. This suggests that the progenitor star 
of SN 2006JC is a WCO Wolf-Rayet star with a C-rich 
envelope and CSM (Figs.Hi-lli). 

Inspection of all the pre supernova models available 
in'Li mongi fc Chieffil ()2006( ) indicates that only massive 
models, i.e., M^s > 4OM0, fulfill the requirements from 
the IR observation and become WCO stars. Moreover, 
these are the only stars in which the chemical composi- 
tions of the mantle and CSM are dominated mainly by 
C with a smaller amount of O (Fig. [5J;d) . 

In stars with initial masses smaller than Mms '^ 35Mq , 
the mass of the He convective core increases or remains 
constant during the core He burning phase. At core 
He exhaustion, a sharp discontinuity of He abundance 
is produced at the outer edge of the CO core. Then, 
the CO core begins to contract to ignite the next nuclear 
fuel while He burning shifts to a shell inducing a forma- 
tion of a convective zone. The He convective shell forms 
beyond the He discontinuity at the outer edge of the 
CO core. Hence its chemical composition is dominated 
by He [X(He)> 0.9]. Because of the short lifetime of 
the advanced burning stages, only a small amount of He 
is burned inside the shell before the presupernova stage 
(Figs.[3^b). Such a behavior is typical for stars in which 
the He core mass remains r oughly constant during core 
He burning (e.g. jNomoto fc Hashimoto 1988). 

In stars with initial masses greater than Afms '^ SSMq , 
on the contrary, the mass loss is efficient enough (10"^ — 



TABLE 4 

Basic evolutionary properties of the 

proc;enitor star. 



Key quantities 


H Burning 


tH [Myr] 


4.64 


Mcc [Mq 


25.80 


Mtot [Mq 


35.40 


to [Myr] 


4.16 


Mho [Mq] 


10.01 


He Burning 




<Hc [Myr 


0.46 


Mho.CC Mq] 


12.56 


Mtot [Mq] 


7.04 


Menv [M0] 


18.80 


X(12C)ce„ 


0.28 


tred [Myr] 


0.07 


twNL Myr (X(He)con) 


0.11 (0.77) 


*WNE Myr (X(He)con) 


0.054 (0.44) 


<wco Myr (X(Hc)ccn) 


0.21 (0.31) 


Advanced Burnings 




Aiexp [yr] 


1.25(+4) 


MHe(max) [Mq] 


16.52 


Mcc (max) [Mq] 


4.83 


Mpe.preSN [Mq] 


1.50 


MpreSN [Mq] 


6.88 


-RproSN [cm] 


3.08(+10) 


M.(He,hell) (Int.-Ext.) [Mq] 


5.262-6.648 


M.(C,hell) (Int.-Ext.) [Mq] 


2.736-4.097 


iwNL(tot) yr 


1.10(+5) 


iwNE(tot) yr 


5.43(+4) 


twco(tot) yr 


2.21(+5) 


twR(tot) [yr] 


3.86(+5) 



10~^Mq yr~^) to uncover the He core and they re- 
duce progressively their mass during the core He burning 
phase. The star enters the WNE Wolf-Rayet stage and 
its subsequent evolution is governed by the actual size of 
the He core. In particular, as the He core progressively 
reduces due to the mass loss, the star tends to behave 
as an initially-lower mass star, i.e., essentially reduces 
its central temperature. This induces the He convective 
core to shrink progressively in mass as well, leaving a 
layer with a variable chemical composition that reflects 
the central abundances at various stages during core He 
burning. When the stellar mass is reduced below the 
maximum extension of the He convective core, the prod- 
ucts of core He burning appear on the surface and the 
star becomes a WCO Wolf-Rayet star. At core He ex- 
haustion. He burning shifts to a shell inducing the forma- 
tion of the convective shell. The convective shell forms 
in the region with variable chemical composition. As a 
consequence, at variance with what happens in stars with 
M^s %j 35Mq, in these stars, the chemical composition 
of the convective shell becomes a mixture of the central 
He burning products. Hence it is mainly composed of C, 
OandHe(Figs.[3td). 

Since all the models above 40 Mq have a similar pre- 
supernova structure, we selected a 40 Mq star as repre- 
sentative of a typical star becoming a WCO Wolf-Rayet 
star. The mass at the presupernova stage (Mpj-cSN) is 
-MproSN = 6.9Mo because of the strong mass loss. We 




Fig. 4. — Presupernova evolutionary path of the progenitor star 
with Mnis = 40AfQ in the Hertzsprung-Russell diagram. 
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Fig. 5. — Temperature and density structures of the presupernova 
progenitor star with Mms = 40Mq . 

underline that the Mms-MproSN relation is highly uncer- 
tain because it strongly depends on many details of the 
stellar evolution (e.g., the mass loss, ov ershooting, rota- 
tion a nd metalhcitv. seeJLa ngcr 198 91: iNugis fc Lamerd 
2000"; 'Mevnct & Maeder'~'2G03; Nomoto et alj |2006|; 
Limongi & Chicffi 2006; Eldridgc & Vink 2006). For this 
reason, for the purpose of this study, we mainly focus on 
a WCO progenitor with MpreSN ~ 6.9 Mq, without pay- 
ing much emphasis on Mms- 

A detailed discussion of the presupernova evolution 
during all the nuclear burning stages is beyond the pur- 
pose of this paper. Hence we report here in Table |4] some 
key properties during the H, He, and advanced burning 
stages. In particular, for the H burning stage we report 
the following quantities: the H burning lifetime (tn), the 
maximum extension of the convective core (Mcc)j the 
total mass (Mtot) at core H exhaustion, the time spent 
as an 0-type star (io) and the He core mass (Mhc) at 
H exhaustion. Here, we assume that the temperature 
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Fig. 6. — Density structures at 100 s after the explosions for 
the models with £51 = 1 (solid), B51 = 5 (dashed), E^i = 10 
(dash- dotted), and £51 = 20 (dotted). 

of the 0-type stars is 33000 K < Teff < 50000 K. For 
the He burning phase we report the following quanti- 
ties: the He burning lifetime (ino), the maximum size 
of the He convective core (Mhc,cc), Aftot at core He 
exhaustion, the maximum depth of the convective en- 
velope (Menv), the central ^^C mass fraction at core He 
exhaustion [X(-'^^C)cen]j the time spent at the red side 
(log Toff < 3.8) of the HR diagram (i^od), and the WNL, 
WNE, and WCO lifetimes (twNL, ^wne, and ^wco, re- 
spectively) - in parenthesis the central He mass fraction 
[X(He)con] when the star enters the WNL, WNE and 
WCO phases. For the advanced burning stage we report 
the following key quantities: the time until the explosion 
(Atcxp), the maximum size of the He core [AfHc(niax)], 
the maximum size of the CO core [Mco(iiiax)], the 
masses of the iron core (MFe,preSN) and the star (MprcSN) 
and the radius of the star (i?preSN) at the presuper- 
nova stage, the final extension in mass of the He con- 
vective shell [Mr(Heshoii)] and of the convective C shell 
[Mr(Cshoii)], and the total lifetimes during the WNL 
[iwNL(tot)], WNE [iwNE(tot)], WCO [%co(tot)], and 
WR [iwR(tot), where iwR = twNL + ^wne + %co] 
phases. 

Figures [4] and [5] show the evolutionary path in the HR 
diagram and the temperature and density profiles at the 
presupernova stage. 

4. HYDRODYNAMICS AND NUCLEOSYNTHESIS 

The SN explosion and explosive nucleosynthesis are 
calculated for the progenitor star with Mpj-eSN = 
6.9Mq. We apply various explosion energies {E51 = 
E/10^^ ergs = 1, 5, 10, and 20) for t he SN explosion 
calculations (e.g., iNomoto et al.l 120061 : iTominaga et al.l 
|20O7b). The hydrodynamical calculation is per- 
formed by means of a spherical Lagrangian hydro- 
dynami cs code with a pie c ewise parabolic method 
(PPM, IColella fc Woodward! \19S4 ) including nuclear 
energy production from the a-network. The equa- 
tion of s tate takes account of the gas, radiation, e~- 
e+ pair (jSugimoto fc Nomotdll975h . Coulomb interac- 
tions between ions and electrons, and phase tran sition 
(|Nomoto et al.lll982t iNomoto fc Hashimotolll988f ). Af- 
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Fig. 7. — Comparison between an escape velocity (dashed) and 
an ejecta velocity for the model with E^i = 1 (solid). 

ter the hydrodynamical calculations, nucleosynthesis is 
calculated as a post-processing with a reaction network 
that includes 280 isotop es up to ™Br (see Table 1 in 
lUmeda fc Nomotol[200l . 

Since the explosion mechanism of a core-collapse SN 
for a massive star with an iron core is still an un- 
solved problem (e.g.. iJanka et al.ll2007[ ). we initiate the 
SN explosion as a thermal bomb. Although there are 
various ways to simulate t he explosion (e.g., a kinetic 
piston, Woos lev fc Weave"rl [l995 ) , it is suggested that 
the explosive nucleosynthesis does not depend sensi- 
tively on the way how th e explosion energy is deposited 
(jAufderheide et all Il991h . We set an inner reflective 
boundary at Mr = IMq and r = 1000 km within the 
iron core and elevate temperatures at the inner bound- 
ary. 

In the spherical symmetry case, for any given progeni- 
tor model, hydrodynamics and nucleosynthesis are deter- 
mined by the explosion energy. During the SN explosion, 
a shock propagates outward inducing local compression 
and heating, triggering explosive nucleosynthesis. Be- 
hind the shock front the matter is accelerated and starts 
moving outward. However, if the progenitor has a deep 
gravitational potential and the explosion energy is low, 
the inner layers begin to fall back due to the gravitational 
attraction. A more compact star and a lower explosion 
energy leads a larger amount of fallback. The fallback 
has a deep implication on the SN nucleosynthesis be- 
cause it decreases the matter ejection, especially, of the 
inner core (e.g., ^^Ni). 

Figure [6] shows density structures at 100 s after the 
explosions when homologously expanding structures are 
reached {v oc r). We find that the fallback takes place 
for the model with E^i = 1 but not for the models with 
E51 = 5, 10, and 20. Figure [7] shows a comparison be- 
tween the escape velocity and the ejecta velocity for the 
model with E^i = 1 and demonstrates that the matter 
below Mr — Mfaii = 3.8Mq will fall back. On the other 
hand, in the models with E51 = 5, 10, and 20, the matter 
above the inner boundary will be ejected. 

The abundance distributions after the explosions are 
shown in Figures [H^-[Hll. In every model, ^^Ni is synthe- 
sized in the innermost layer {Mr < A/pcprcSN = 1.5A/q) 
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Fig. 8. — Abundance distributions after the explosions of the progenitor star with Mms = 4OM0. The explosion energies are (a) E^\ = 1, 
(b) £51 = 5, (c) E51 = 10, and (d) E51 = 20. 



TABLE 5 
Nucleosynthesis properties of the explosion 

MODELS WITH Mp^eSN = 6.9M0 (Mms = 40Mq). 



Explosion energy 

[10^1 ergs] 






Msi 



M(56Ni) 



Mcut 



1 

5 
10 
20 



1.8 
2.1 
2.3 
2.5 



2.1 
2.5 
2.7 
3.0 



0.5 
0.6 
0.7 



1.8 
2.0 
2.3 



due to the low electron fraction. Thus, we can estimate 
the maximum amounts of synthesized ^^Ni for given en- 
ergies. The ^^Ni-rich layer extending to M^. = Mpo 
[where ^(^'^Ni) = ^(^^Si)] and ^^Si-rich layer extending 
to Mr = Afsi [where Xi^^'ax) = X{^^Q))\ expand farther 
in the models with higher E because the temperature 
achieved is higher in the outer layer for higher E. Afpc 
and Mgi for each model are summarized in Table [5l 

^•^Ni is synthesized at M^ < Msi- Since Msi < Mfaii in 
the model with £'51 = 1, the model is likely not to eject 
^^Ni. On the other hand, the models with E51 = 5, 10, 
and 20 can eject all synthesized ^^Ni because the fallback 
does not occur. The total amounts of synthesized ^^Ni 
for the models with E^i ~ 5, 10, and 20 are summarized 
in Table [3 



5. LIGHT CURVE 

The energy source of the LC of SN 2006jc is still un- 
der debate. The possible sources include the '''^Ni-^^Co 
decay like Type I SNe and the ejecta-CSM interaction 
like Type Iln SNe. However, both scenarios have the fol- 
lowing problems. In the case of the ^^Ni-^^Co decay, the 
7-ray photon and positron emitted from the ^^Ni-^^Co 
decay are absorbed by the SN ejecta and the absorbed 
energy is thermalized. Thus, the spectra would show 
a blackbody-like continuum as normal Type I SNe do. 
However, the spectra of SN 2006jc do not resemble those 
of normal Type I SNe but sh ow a bri ght blue c ontinuum 
in early epochs (jFolev et al.1 '2007: S mith et al. ! 2008). In 
the case of the ejecta-CSM interaction, the kinetic energy 
is transformed to an X-ray emission via bremsstrahlung 
radiation, and then converted to UV, optical, and IR 
emissions. Thus, it is difficult to explain that the X-ray 
luminosity is much fainter than the optical luminosity 
unless the optical depth for the X-ray emission is much 
higher than that for the optical emission. Another prob- 
lem with the ejecta-CSM interaction model is that the 
X-ray LC is not synchronized with the bolometric LC. 
In addition, the LC powered by the ejecta-CSM inter- 
action usually has a long-term plateau (e.g., SN 1997cy, 
iTuratto et al.l [2000 1. Thus, we assume that the LC is 
powered by the ^^Ni-^^Co decay. 



5.1. Radioactive Decay Models 

The bolomctric LC of SN 2006jc is constructed from 
the UV, optical, and IR observations as described in 
§ [2] The estimated peak bolomctric magnitude of 
SN 20 06.ic is M = -18 .4, being as bright as SN 2006aj 
(e.g., iPian et all l2006f) . Thus, it is speculated that 
the ejected amount of ^^Ni [M( ^^Ni)] is similar to 
SN 2006ai, i.e. , Mi^'^Ni) ~ 0.2Mq (IMazzah et al.ll2006L 
I2007t iMaed'aet al. 2007) . According to § H the mod- 
els with E51 = 5, 10, and 20 can eject a large enough 
amount of ^^Ni, while the ^^Ni production of the model 
with E51 = 1 is too small. 

The spherical explosion models with E51 — 5, 10, 
and 20 yield too much M(^^Ni) because of no fallback. 
However, no fallback is a consequence of the assump- 
tion of the spherical symmetry. The fallback takes place 
in an aspherical explosion even with a high explosion 
energy and thus the aspherical explosion may well de- 
crease M(^^Ni) and increase the central remnant mass 
Mr.^ (iMaeda fc Nomotol 120031: iTominaga eFail l2007al : 
iTominagal l2007f K Therefore, assuming that aspheri- 
cal fallback takes place in the high-energy models with 
£■51 = 5, 10, and 20, we estimate the amount of fallback 
to yield M(^^Ni) ~ O.2M0 and then the ejected masses 
for the models as Moj — MpreSN — -^rem- As a result, the 
sets of Mrcm, M^j and E are derived to be {M^b^/Mq, 
Mej/Mg, E51) = (1.8, 5.1, 5), (2.0, 4.9, 10), and (2.3, 
4.6, 20). 

Applying the homologous density structures of the 
models (Fig. [S]), we synthesize bolomctric LCs for the 
models with E^i = 5, 10, and 20 using the LTE radiation 
hydrodynamics code and the gray 7-ray transfer code 
(jlwamoto et al.ll2000f ). In the radiative transfer calcula- 
tion, the electron scattering is calculated for the ioniza- 
tion states solved by the saha equation and the Rosseland 
mean opacity is approximated wit h an empirical rela tion 
to the electron-scattering opacity (|Deng et al.ll2005f ). 

The peak width (r) of the SN LC depends on 
the ejected mass Mej, explosion energy E, opacity 
K, density structure, and ^^Ni distribution, as r oc 
Ak^/'^M^/'^E-^/'^ ()Arnettlll98l ). where A represents the 
effects of the density structure and the ^^Ni distribution. 
Here, we assume for sake of simplicity a uniform mixing 
of ^^Ni in the SN ejecta. Also, the density structures 
after the SN explosions with various E are analogous. 
Thus, the dependence on A is negligible and we investi- 
gate the LC properties depending on k, M^j and E. The 
synthetic LCs obtained for the models with (Mej/M©, 
E51) = (5.1, 5), (4.9, 10), and (4.6, 20) are shown in Fig- 
ure [51 Figure [H] also shows the multicolor and bolomctric 
LCs of SN 2006JC. 

5.2. Comparison with Observations 

The period of SN 2006jc is divided into four epochs 
depending on the available observations: (1) UV and op- 
tical photometries at i < 50 days, (2) optical and NIR 
photometries at i ~ 50 — 80 days, (3) optical photometry 
at t - 80 - 120 days, and (4) optical, NIR, and MIR 
photometries and NIR spectroscopy at i > 120 days. 

(1) At i < 50 days, the IR contributions to the bolo- 
mctric luminosity may well be small because the IR con- 
tribution is only ~ 10% at i ~ 50 days. Thus, the peak 
bolomctric luminosity derived from the UV and optical 
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Fig. 9. — Comparison between the synthetic LCs for the mod- 
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5.1Mq (green line), B51 
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20 and M^ 



4.6Mq [red li ne) and the LCs of SN 2006jc (L^jy -|- Lppt : cyan 
filled circles , llmmler et aP l2008l : IKawabata et al.l I2007al . l2008l : 
' Anupama e t al. 2008; J^m.c stf'^ < 3 X lO ^^Hz): red filled circles, 
Arkharov et al. 2006; Di Ca rlo et al.ll20071; L^^ol^■ b lue filled circles, 
Immlcr et al. 2008; Kawabata et al. 2007a. |2008|: lAnupama et al.l 
2008; Arkharov ct al. 2006; Di Carlo et al.l 120071 ; Lra h„t(i^ < 
3 X lO^'^Hz): magenta filled circle, Sakon ct al. 200^; Li^{u < 
3 X lO^^Hz): magenta open circle, Sakon ct al. 2007). The lumi- 
nosities denoted by the circles without errorbars are considerably 
contributed by the UV luminosity estimated as Luv = 0.8Lopt. 

fluxes is reliable (§ 12. ip . If the bolomctric LC peaked at 
the discovery, the peak luminosity is reproduced by the 
56Ni_56co decay of M(^6j^i) = 0.2271//©. The rapid de- 
cline after the peak prefers such high-energy models as 
(Mej /Mo, E51) = (4.9, 10) and (4.6, 20). 

(2) At t ~ 50 - 80 days, the IR contribution to the 
bolomctric luminosity increases from ^ 10% at i = 49 
days to ^ 70% at i = 79 days. The contribution of 
iiR,cst(i/ < 1.3 X lOi^Hz) to Lin,est{iy < 3 x lO^^Hz) 
changes from ~ 30% at i = 49 days to ^ 40% at t = 
79 days. At this epoch, the optical and IR emissions 
contribute to the bolomctric luminosities. Combining the 
IR brightening and the optical decline, the bolomctric LC 
including Lyv {— 0.8Lopt) declines slowly. Such a slow 
decline is consistent with the models of {Mcj/Mq, E^i) 
= (4.9, 10) and (4.6, 20). 

(3) At i ~ 80 - 120 days, NIR photometries are not 
available. The decline of the optical luminosity at this 
epoch is more rapid than at t < 80 days. Such a rapid 
decline of the optical LC can not be reproduced by the 
56]\jj_56(^Q decay. However, the bolomctric LC may well 
decline more slowly than the optical LC because the IR 
emission dominates in the bolomctric luminosity. 

(4) At < > 120 days, NIR photometries are available 
continuously and optical p hotometries are availa ble at 
t = 120, 140, and 195 days (jKawabata et al.ll2008^ . The 
contribution of the optical emission to the bolomctric lu- 
minosities is negligible (^ 3% at t = 120, 140, and 195 
days). At this epoch, the contribution of iiR,est(^ < 
1.3 X 10^'^Hz) to the total luminosity increases from 
~ 60% at t = 127 days to - 80% at t = 215 days. The 
estimated IR luminosity is consistent with the luminosity 
emitted from the hot carbon dust at i = 215 days (S I2.2p . 
Since the dust temperature decreases with time, the ratio 
of the MIR luminosities to the NIR luminosities becomes 
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Fig. 10. — Comparison between the synthetic X-ray LCs con- 
tributed from the total emission {solid line), the SN ejecta {dashed 
line), a nd the CSM (dotted line) and the X-ray LC of SN 2006jc 
{circles, llmmler et aL|[200a) . 

larger with time. Therefore, the amorphous carbon emis- 
sion model reasonably estimates the IR luminosity due 
to the hot carbon dust at t < 215 days. The model with 

Afoj = 4.9Mq, E51 = 10, and M{^^m) = 0.22Mq re- 
produces well the LC decline at i > 120 days and the 
IR luminosities due to the hot carbon dust at t = 215 
days^^. Therefore, we conclude that the hypernova-like 
SN explosion model with Mcj = 4.9Mq, E51 = 10, and 
M(^^Ni) = O.22M0 is the most preferable model among 
the exploded models of a WCO Wolf-Rayet star with 

MproSN = 6.9Mq. 

6. INTERACTION WITH CIRCUMSTELLAR MATTER 

X-rays from SN 2 006.ic were detected by the Swift and 
Chandra satellites (jlmmler et al.ll2008D . The X-ray de- 
tection indicates that the expanding SN ejecta collides 
with the CSM. 

We calculate X-ray emission from the ejecta-CSM in- 
teraction for the SN model with {M^j/Mq, E51) = 
(4.9, 10), and estimate the CSM density structure on 
the b asis of a comparison wi th the observed X-ray LC 
(e.g., ISuzuki fc Nomotolll995D . The observed X-ray lu- 
minosities es t imate d with the distance of 24 Mpc in 
llmmler et all ()2008f ) are scaled using 25.8 Mpc. 

We adopt a CSM density profile characterized by a 
power-law of p = Po(^/''o) " and assume that the inter- 
action starts at a distance r = 3 x 10^° cm. The parame- 
ters po, ro, and n are determined so that the ejecta-CSM 
interaction reproduces the observed X-ray LC. 

The interaction generates reverse and forward shock 
waves in the SN ejecta and CSM, respectively. Both re- 

^® After submission of this paper, the IR observation at t = 425 
days was presented by IMattila et all II2008I) . They estimated the 
total luminosity and the dust temperature as L ~ 1.2x10'*'' ergs s~* 
and T = 520 K, respectively. Our model with M^j = 4.9Mq, 
E51 = 10, and M(56Ni) = O.22M0 predicted the luminosity of 



^bol ~ 3 X 10'^'^ ergs s~^and the dust temperature of T ~ 200 K 
IjNozawa et al.l [20081 ) at i ~ 430 days, which are lower than the 
observations. This suggests that the observed IR emissions at t = 
425 days may originate not only from the newly-formed dust in the 
SN ejecta heated by the ^^Ni-^^Co decay but also from the light 
echo of the CSM dust. 
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Fig. 11. — (a) Density, (b) velocity, and (c) temperature struc- 
tures of the SN ejecta and the CSM at t = day (red), 10 days 
{green), 20 days {blue), 50 days {magenta), 100 days (cyan), and 
200 days {black). The coordinate is the Lagrangian mass with the 
contact discontinuity between the ejecta (left) and the CSM {right). 



gions are heated by the shock waves and emit X-rays. In 
such a compact star, because the density in the shocked 
SN ejecta is higher than that in the shocked CSM, the 
emitted X-rays from the shocked SN ejecta are more lu- 
minous than those from the shocked CSM (Fig. ITU)). 

Figure (TOj shows the synthesized X-ray LC for po = 
2.75 X 10"^^ g cm"3 and n = for r < 2.2 x 10^*^ cm and 
n — 6 for r > 2.2 x 10^^ cm, i.e., for a flat (inside) and 
steep (outside) CSM density profile of p = 2.75 x 10"^^ 
g cm~3 foj. r <2.2x lO^^ cm and 2.75 x 10-^^{r/2.2 x 
lOi'^cm)-'^ g cm-3 foj. r > 2.2 x 10^^ cm. The total mass 
of the CSM is 1.2 x IO^'^Mq to reproduce the peak of 
the observed X-ray LC and the subsequent decline. 

The density, velocity, and temperature structures and 
their evolutions are shown in Figures ITTk bc. The velocity 
of the reverse shock is w '^ 3.8 x 10^ km s~^. The reverse 
shock reaches ~ 5.3 x IO^^Mq from the outer edge of 
the SN ejecta at t = 200 days and heats up the swept-up 
SN ejecta. The temperature behind the reverse shock is 
higher than 10^ K where dust cannot newl y form and the 
dust formed in the SN ejecta is destroyed (jNozawa et al.1 
l2008f ). Our calculation does not show the formation of a 
cooling shell. This is because the CSM interaction is so 
weak to emit X-ray of '^ 3 x 10^^ ergs s~^. If the bolo- 
metric luminosity is powered by the CSM interaction, 
i.e., if the CSM interaction emits as high luminosity as 
~ 8 X 10'*^ ergs s~^, the cooling shell might form and 
thus the dust formation might be possible. Further de- 
tailed studies, however, are required to confirm the dust 
formation behind the reverse shock. 

Such a flat density profile of the inner CSM implies 
that the stellar wind was not steady because the steady 
wind should form the CSM of p oc r~^. This circumstel- 
lar environment might have been formed by a variable 
mass- loss rate M and/or a variable wind velocity v^j. 
For example, assuming that the stellar wind blew with 
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a constant v^, = 3,500 km s^^ for two years before the 
explosion, the mass-loss rate must have changed from 
1 X 10~^ to 2 X 10~^^ Mq yr~^ in two years. Such a 
drastic change of the mass-loss rate and/or the wind ve- 
locity is consistent with the fact that the progenitor of 
SN 2006JC was surrounded by the matter ejected by the 
LBV-like event two years before the explosion. 

7. CONCLUSIONS 

We present a theoretical model for SN 2006jc whose 
properties are summarized as follows. 

(1) WCO progenitor and Dust Formation: The 

progenitor is a WCO Wolf-Rayct star whose total mass 
has been reduced from M^s — 4:0Mq to as small as 
MproSN = 6.9AfQ. The WCO star model has a thick 
C-rich envelope and CSM. This is consistent with the for- 
mation of amorphous carbon grains in the S N eject a and 
the C S M suggested by AKA fil observations (jSakon et al.l 
l2007t ). iNozaw-a et~alT(j2008[ ) have calculated dust forma- 
tion in the WCO star explosion model and shown that 
carbon dust is formed in the C-rich layer at t ^ 50 
days. This is much earlier than the dust formation after 
i ~ 1 yr in Type 11 SNe (SNe II), because of the much 
smaller Mej in the WCO star than SNe II. According to 
the models in iLimongi fc Chieffil ()2006l ). the stars with 
Myns > 4:0Mq typically become WCO stars to form a 
thick C-rich layer and CSM. This limiting mass, however, 
is still uncertain and strongly depends on many details 
of the stellar evolution. The early dust formation in the 
SN ejecta and the CSM suggests that the progenitor of 
SN 2006JC is a massive star becoming a WCO Wolf-Rayet 
star. 

(2) Explosion and Bolometric Light Curve: The 
multicolor LCs of SN 2006jc show peculiar evolutions, 
e.g., a rapid decline of the optical LC and brightening 
of the IR LC. These can be interpreted as an ongoing 
dust formation. Assuming the absorbed optical light is 
re-emitted in the IR band, the bolometric LC is con- 
structed as a summation of Lyv (= O.SLopt), ^opt, and 
LiR. By calculating the hydrodynamics, nucleosynthe- 
sis, and the bolometric LC for the SN explosion with 
the various explosion energies, -E51 = 1,5, 10, and 20, 
wc find that the hypernova-like SN explosion model with 
Afej = 4.9M0, E51 = 10, and A-f (56Ni) = O.22M0 best 
reproduces the bolometric LC of SN 2006jc with the ra- 
dioactive decays. Also, the temperature evolution of the 
carbon dust heated by the ^^Ni-^^Co decay reasonably 
well explains the IR observations for i ^ 50 — 220 days 
(jNozawa et al.ll2008[ ). 

(3) CSM Interaction and X-ray Light Curve: Ap- 
plying the model with A/cj — 4.9A/0 and £'51 = 10, we 
calculate the ejecta-CSM interaction and the resultant 
X-ray LC. We derive the CSM density structure to re- 
produce the X-ray LC of SN 2006jc as p = 2.75 x lO'^^ 
g cm-3 for r < 2.2 X 10^^ cm and 2.75 x 10-^^{r/2.2 x 
10i^cm)-6 g cm-3 foj. ^ > 2.2 x 10^^ cm. The flat density 
distribution in the inner CSM indicates a drastic change 
of the mass-loss rate and/or the wind velocity that is 
consistent with the LBV-like event two years before the 
explosion. 



8. DISCUSSION 

LBV connection: Our model does not take into ac- 
count the LBV-like event that occurred two years before 
the explosion. The first reason is that the mechanism 
of the outburst is still unclear. The second reason is 
that, at least in the framework of the current under- 
standing of standard stellar evolution, the envelope of a 
massive star practically freezes out after core He exhaus- 
tion (i.e., about 10000 years before the explosion) due to 
the more rapid evolution of the core than the envelope. 
In addition, it is interesting to note that, and this is a 
confirmation of the theoretical expectation, there is no 
observational evidence that any Wolf-Ray et star has ever 
unde rgone such a luminous outburst ( Humphrevs et aLl 
Il999t ). Hence, there is no specific reason to associate the 
occurrence of a LBV-like outburst to the presupernova 
evolution. Future studies on the mechanism of the out- 
burst are required to firmly conclude the origin of the 
LBV-like outburst. It would be possible that a possi- 
ble binary companion star could undergo the LBV-like 
outburst. 

Fallback: According to our hydrodynamics and nu- 
cleosynthesis calculations, in the spherically symmetric 
models with E51 > 5, the fallback does not take place and 
thus the amount of synthesized ^^Ni is much larger than 
Af (^^Ni) = O.22Af0 which is required to power the LC of 
SN 2006JC ('i) l5.1|) . In the aspherical explosions, however, 
the fallback takes place even for E51 > 5. In this paper, 
we assume the fallback even for the models with E51 > 5 
and derive the amount of fallback to yield the appropri- 
ate amount of ^^Ni. To justify the above assumption, we 
calculate an aspherical explosion induced by a jet with an 
opening angle of9 = 45° and an energy deposition rate of 
E = 3 X 10^2 ergs s"^ (jTominaga et al.ll200 7atlTominagal 
|2007| ). The jet-induced model realizes an explosion with 
Afoj ~ 4.9Af0, E51 ~ 10, and Af (S^Ni) - O.22A//0 that is 
consistent with the adopted model. We note that an as- 
pherical radiative transfer calculation is required to con- 
firm that the jet-induced explosion model can reproduce 
the LC of SN 2006JC. 

Light curve models: The model with Mcj — 4.9A'f0, 
E51 = 10, and M{^'^m) = O.22Af0 is not an unique 
model to reproduce the bolometric LC of SN 2006jc. In 
the case of usual SNe, the velocities of the absorption 
lines can disentangle the degeneracy of Afej and E by 
means of the compariso n with the photospheric veloci- 
ties (e.g.. pTanaka et aLl 2008). However, the spectra of 
SN 2006JC are dominated by He emission lines and the na- 
ture is unclear. Thus we cannot fully resolve the degener- 
acy. Since the LC shape is proportional to M-' E^^^'^, 
the model with a larger Mcj requires a higher E. An 
explosion of the progenitor star with a larger AfpioSN 
reproduces the LC of SN 2006jc with a higher E and 
the X-ray LC with a lower CSM density. On the other 
hand, an explosion of the progenitor star with a smaller 
AfproSN reproduces the LC with a lower E, e.g., an ex- 
plosion with Afoj = 1.5Af0 and Er,i = 1 can explain the 
LC shape of SN 2006jc. However, such low E explosions 
suppress explosive nucleosynthesis and enhance the fall- 
back. As a result, the ^^Ni production is reduced for 
a small A^prcSN- If the LC of SN 2006jc is powered by 
the ^^Ni-^^Co decay, the bright peak indicates a larger 
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amount of ^^Ni production {^ 0.22M(7,) than a normal 
SN[M (56Ni) ~ O.O7M0, e.g., SN 1987A. lBlh^ikov et all 
120001 1 . Therefore, SN 2006jc is likely a more energetic ex- 
plosion than a normal SN with E51 ^ 1. 

Dust formation: We assume that the energy source 
of the LC of SN 2006jc is the ^^Ni-Seco decay. This 
consistently explains the formation of carbon dust at 
the early epoch {t ^ SOdays) and the dust tempera- 
ture at t '^ 200days (jNozawa et al.l 12008') . In this sce- 
nario, however, the origin of the bright blue continuum 
remains an unsolved problem (e.g.JPastorello et ani2007l : 
ISmith et al.]l2008t llmmler et al.l l2008[). Such a spectrum 
might be explained by the ejecta-CSM interaction. In 
this scenario, however, the fine tunings are required to 
reproduce the bolometric LC; most of the X-rays are ab- 
sorbed and converted to the optical luminosity, which 
only a small fraction of the X-rays are emitted with 
changing the fraction from 10~^ at i ~ 30 days to 0.1 
at i ~ 180 days. Moreover, the formation of carbon dust 
with two temperatures would not be explained. Since 
both scenarios are inconclusive so far, further investiga- 
tions may give important implications on the emission 



mechanism of SN 2006jc. 
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